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ABSTRACT 

We prove that a class of solutions to Einstein's equations — originally discovered by G. 
C. McVittie in 1933 — includes regular black holes embedded in Friedman-Robertson- Walker 
cosmologies. If the cosmology is dominated at late times by a positive cosmological constant, 
the metric is regular everywhere on and outside the black hole horizon and away from the 
big bang singularity, and the solutions asymptote in the future and near the horizon to the 
Schwarzschild-de Sitter geometry. For solutions without a positive cosmological constant the 
would-be horizon is a weak null singularity. 
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1 Introduction 



Finding solutions to Einstein's equations describing anything beyond the simplest and most 
symmetric configurations of matter or gravity is a hard and uncertain affair. When available 
they deserve our full attention, as they can provide us with new insights into the nature of 
gravity beyond the linearized regime or be of use in describing some objects of astrophysical 
relevance. For both reasons any solution potentially describing a black hole embedded in an 
expanding universe is of considerable interest. 

It is thus a bit surprising that a proper understanding of a class of solutions found over 
70 years ago by McVittie [1] is still lacking. These solutions have many of the features one 
would expect of a black hole embedded in a Friedmann- Robertson- Walker (FRW) cosmology: 
they are spherically symmetric, parametrized by a function a{t) and a mass parameter m, 
reduce to FRW cosmology with scale factor a(t) at large radius, and reduce to known black 
hole metrics or standard FRW cosmology in all the appropriate parametric limits. These 
properties make them not only interesting in their own right as non-linear solutions, but 
potentially significant and physically relevant for describing real gravitating objects or holes 
in our universe. On the other hand there is no accretion (the mass m is constant) — an odd 
property for a physical black hole in a universe full of matter or radiation. 

The most recent and thorough examination of the properties of the McVittie metric was 
attempted in a series of three papers by Nolan [2] , which review past work on the metric and 
describe many of its features. Because it is spherically symmetric and asymptotically FRW, 
the McVittie metric can be used to describe the external fields of finite size objects [2, 3, 4], 
or exteriors of bubbles separating different regions of spacetime [5, 6]. In these applications 
the McVittie metric is replaced by a different geometry at small radius. But — as Nolan 
reviews at length — ^the literature on the McVittie solution proper is riddled with basic errors 
and contradictory statements, which has probably contributed significantly to its relative 
obscurity. 

Unfortunately Nolan himself erred on some points that are critical for the correct in- 
terpretation of the solution. Specifically, he argued (but didn't prove or derive) that the 
would-be null black hole horizon of the McVittie metric is at infinite distance and therefore 
constitutes a null boundary rather than a horizon, and hence that the metric outside this 
surface is geodesically complete and cannot describe a black hole at all. Much of the rest of 
his analysis relied on this assertion, which as we will demonstrate explicitly is incorrect. The 
null surface is at a finite distance and therefore renders the standard form of the McVittie 
metric geodesically incomplete, a conclusion that entirely changes the interpretation given 
in [2] of the global structure of the spacetime and vahdates the black hole interpretation in 
at least some cases. 

The cases for which the black hole interpretation is valid occur when the McVittie scale 
factor asymptotes to de Sitter space^ (a(t) e^°*). The late-time exponential expansion 
turns out to dilute the FRW fluid(s) near the horizon rapidly enough that all curvature 

^If one takes a negative cosmological constant the cosmology must asymptote to an open FRW universe 
which originates from a big bang singularity and terminates in a big crunch singularity (see, e.g. [7] for a 
discussion). Similarly, in collapsing FRW universes with any value of the cosmological constant there would 
be spacelike or null singularities in the future. We will not consider these cases further. 
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invariants reduce to their corresponding values for a pure de Sitter-Schwarzschild hole. While 
we are unable to explicitly construct the analytic extension across the horizon, using the null 
ingoing geodesies we can show that the metric is nondegenerate and analytic near the horizon, 
which implies the existence of such an extension. In contrast, for generic power law scale 
factors (a(t) ~ at late times) the expansion is not fast enough to sufficiently dilute the 
stress-energy on the would-be horizon, rendering it weakly singular. We explicitly construct 
a higher derivative curvature invariant that diverges on the null surface in this case, implying 
the existence of infinite tidal forces there (and potentially large quantum corrections arising 
from higher curvature invariants in the effective action). 

To conclude this introduction we list the sahent features of a spatially fiat McVittie 
metric characterized by the scale factor a{t) of the asymptotic FRW region and a mass m: 

• for m = 0, reduces to standard homogeneous and isotropic FRW cosmology with scale 
factor a{t) and zero spatial curvature; 

• for H{t) = dt\na(t) = constant, reduces to a Schwarzschild or de Sitter-Schwarzschild 
black hole of mass m and Hubble constant H; 

• has a spacelike and inhomogeneous big bang singularity in the past (for generic ex- 
panding a{t)); 

o has a null or spacelike FRW future infinity at large radial distance r and late time t; 

o has a surface r = r-,t = oo which is null and at finite affine distance along ingoing 
null geodesies from any point in the bulk; 

o this null surface is a soft, null naked singularity in an FRW spacetime if the FRW 
Hubble constant H(t) = a{t)/a{t) goes to zero at late times; 

• the addition of any positive cosmological constant (so that \iva.t-^^H{t) = Hq > 0) 
eliminates this singularity; 

o all curvature invariants on the null surface in this case are exactly equal to their values 
on the horizon of a Schwarzschild-de Sitter black hole of mass m and Hubble constant 
Ho>0; 

o therefore, at least in the case Hq > the McVittie metric describes a regular (on the 
horizon) black hole embedded in an FRW spacetime. 

The items marked o are new results of the work reported in this paper, correcting existing 
claims in the literature. 

2 McVittie Retrospective 

We begin with a brief review of the McVittie solutions. In their simplest form they have 
zero spatial curvature in the asymptotically FRW region, but can be easily generalized to 
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include non-zero positive or negative spatial curvature [1] or electric charge [8, 9]. We do not 
expect that the spatial curvature of the FRW geometry to significantly alter the behavior 
of the metric near a mass source as long as the gravitational radius of the mass m. or the 
spatial extent of the region occupied by it, whichever is larger, is smaller than the radius of 
curvature.^ Since this is presumably the case for astrophysical masses, we specialize to the 
case of zero spatial curvature. 

The solution is given by the metric 

ds'' = -(^)'^^' + (1 + ^J)'a\t)dx^ , (1) 

where a{t) is the asymptotic cosmological scale factor, 

m . . 

tji/Gm is the mass of the source, and using spatial translations we have chosen x = as the 
center of spherical symmetry. This is an exact solution of the field equations of Einstein's 
General Relativity for an arbitrary mass m provided that a{t) solves the Friedmann equation 

^H{tf = SnGNpit) , (3) 

where p is the energy density T°o and H — a/a. Surprisingly the energy density is constant 
along slices of constant t. It scales with the cosmic scale factor and controls the overall 
expansion rate of the universe exactly as in a standard FRW geometry with scale factor 
a and Hubble parameter H. In fact we can linearize the McVittie solution by considering 
the limit /i ^ 1, in which case the metric (1) reduces to a perturbed FRW universe in 
Newtonian gauge, ds'^ — — (1 — 2^)(ii^ + a^(l + 2-^^)dx^. Starting from this we could 
use linear superposition to build an arbitrary 'multicentered' solution. The result is just a 
general perturbed FRW cosmology. 

The pressure on fixed t-slices of the geometry (1) is not homogeneous. After a slight 
rewriting of the formula in e.g. [11], 

and so has two contributions: a homogeneous term oc H^, and an inhomogeneous term 
oc H. To understand the role of the second term, note that (except in the case of pure 
cosmological constant H = 0) one expects the mass to break the homogeneity of the stress- 
energy on spatial slices. It should pull matter in from the FRW fiuid around it, making the 
energy density inhomogeneous. This does not happen in the McVittie solution: the energy 
density (3) is a function of the cosmic time alone. ^ Therefore something must cancel the 

^In the first of references [2] , Nolan concludes that open McVittie solutions are not good solutions de- 
scribing masses in expanding FRW universes. We suspect that this inference arises from excessively strong 
assumptions. 

^We call this time coordinate the 'cosmic time' since it reduces to the usual comoving FRW time when 
the mass source is absent, and asymptotes to it far away from the source when the mass does not vanish. 
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gravitational attraction of the mass, and that non-gravitational balancing force is provided 
by the gradient of the pressure (4). 

The McVittic solution has a curvature singularity at /i = 1 (which can be seen by 
evaluating e.g. the Ricci scalar R = 12H'^ + Gj^H) where the pressure goes to infinity. As 
we will see this singularity is spacelike, extends all the way to spatial infinity, and should be 
viewed as a cosmological big bang singularity. 

One may wonder what the status of this solution is among candidate metrics describing 
black holes embedded in expanding FRW universes. Nolan [2] offers a theorem proving that 
the McVittie solution is the unique solution describing the field of a spherically symmetric 
mass in a spatially fiat asymptotically FRW cosmology. However, his theorem relies on 
rather strong assumptions, among which are that i) the mass is constant, and that ii) there 
exists a gauge in which energy density is spatially homogeneous despite the gravitational 
field of the central mass. Evidently these conditions can only be attained with the help of 
the non-gravitational forces encoded by the inhomogeneous pressure (4). 

These assumptions are by no means sacrosanct — on the contrary, they are very restrictive 
and not very well-motivated physically. Astrophysical black holes generally accrete matter 
at some non-zero rate. Even solutions with no accretion should exist without the necessity of 
introducing non-gravitational forces, for example by assigning the matter an initial velocity 
distribution such that it fiows radially away from the black hole at or above escape velocity. 
In fact there exist special solutions in the "Swiss cheese" class [13] where a region of fiat 
FRW dust-dominated cosmology is excised and replaced by a void containing a black hole 
with a mass equal to that of the would-be mass of the excised region [12]. One can imagine 
generalizing such solutions further by adding spherical dust shells to the void, each with a 
nonzero peculiar velocity so that e.g. it is marginally bound to the black hole and hence 
never falls in. Presumably one could take a continuum limit of many such shells, finding a 
solution where the central mass remains constant, with zero pressure, but with the energy 
density a function of radius and time. 

Thus, the McVittie solution should be thought of as a special case of a larger class of 
geometries describing masses in FRW: McVittie is the special class where the mass parameter 
is a constant and the energy density is homogeneous, and its inhomogeneous pressure is the 
necessary and sufficient price one pays for these features. 

The initial Big-Bang singularity is absent when H = 0, and in fact the geometries (1) 
reduce to the Schwarzschild or Schwarzschild-de Sitter solutions (for H = and H ^ 
respectively). The hypersurface /i = 1 is perfectly regular in those cases, being the event 
horizon in the Schwarzschild case, and a spacelike surface inside the event horizon in the 
Schwarzschild-de Sitter geometry. Their black hole singularities remain censored by the 
event horizon, and cannot be seen by exterior observes. 

3 Coordinate Covers 

In the case a{t) = 1 the McVittie solution reduces to a black hole in fiat space; setting a = 1 
in the metric (1) gives the Schwarzschild solution in isotropic coordinates. These coordinates 
have the unfortunate feature that the coordinate I a; I covers the exterior of the back hole twice: 
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m/2 < |x| < oo covers the same region — the exterior of the hole — as < \x\ < m/2. 

We will use another coordinate choice which more closely imitates the familiar (static) 
form of the Schwarzschild or Schwarzschild-de Sitter metric. The new radial coordinate is 
defined by [2] 

r = {l + fxfa{t)x, (5) 

where r = |r| turns out to be the 'spherical area' coordinate. This is 'one-half of the 
coordinate transformation which goes from cosmological to static patch coordinates in de 
Sitter space. The analogue of the static patch time is not so easy to find, because in general 
the true time dependence of the McVittic metric makes the integration harder. One can 
show that a time coordinate analogous to static time does exist, but we will not need it for 
our analysis and will use t as a time variable. 

The equation relating a(t)|a;| to r is quadratic, so the coordinate transformation (5) 
actually defines two separate branches: 

To check which of these is physically relevant we note that in the limit r — )■ oo, the product 
a(t)|a;| — )■ on the "+" branch. This branch covers the region of space between a(t)|x| = 0, 
where the curvature blows up at a(t) = for any finite x, and a(t)|x| = m/2, which by (2) 
corresponds to /i = 1 and is another spacelike curvature singularity, the McVittie big bang. 
Since this patch terminates on spacelike curvature singularities both in the past and in the 
future, we will ignore it hereafter and specialize to the "-" branch, where r — ?> oo implies 
\x\ — )■ oo at any finite a{t), so that the geometry asymptotes at large r to an FRW region. 

We are most interested in the solutions which contain both perfect fiuid (or fiuids) and 
positive cosmological constant components. If the fiuid obeys the null energy condition, these 
solutions will be dominated by the cosmological constant at late times. One might expect 
based on the maximally extended Schwarschild-de Sitter geometry that this class of McVittie 
solutions has a maximal extension given by the periodic repetition of the fundamental block 
we will describe (for details, sec e.g. [14]), but we will not pursue this further here. 

If we swap X for r defined in the first of Eq. (5), the McVittie metric (1) becomes 

ds^ = -fdt'' - , dr dt + -J^^ + r^dn2 , (7) 

^Jl - 2m/r 1 - 2m/r 

where / = 1 — 2m /r — H{t)'^r'^. When H = const, this is the Schwarzsczhild-de Sitter metric 
in coordinates which are analogous to outgoing Eddington-Finkelstein coordinates for a flat 
space Schwarzschild black hole (see Figure 1). 

4 Causal Structure 

We will now proceed to use the coordinates (7) to analyze the causal structure of the McVittie 
solution, focusing on the case where p consists of one or more perfect fluid components 
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Figure 1: The section of the Schwarzschild-de Sitter conformal diagram, which is covered by 
our coordinates r,t when H = Hq = const. Black sohd hues are the surfaces of constant t, 
and the green arrow represents an ingoing observer, which crosses these surfaces en route to 
the bold red dashed line r = r_, t = oo, which is the null apparent horizon. This means that 
the coordinate t is a good time coordinate for ingoing observers. 



satisfying the null energy condition plus a positive cosmological constant (so that limf_^o H = 
oo, and limt_>oo H = Hq > 0) ■ The metric function / = 1 — 2m /r — H{t)'^r'^ and its positive 
roots will play important roles in what follows. 

There are five important surfaces that delineate the causal structure. 

• A surface at r = 2m, t = finite. This is the curvature singularity at /i = 1 mentioned 
earlier (c./. e.g. the Ricci scalar R = 12H^ + 6H/^1 — 2m/r). The metric (7) shows 
immediately that it is a spacelike 3-surface: since r is fixed to 2m, we have dr = and 
9u = —f = 4:m'^H'^{t) > 0, and the sphere has finite radius r = 2m. This surface lies 
in the causal past of all spacetime points in the patch of the metric we are studying, 
so we refer to it as the big bang. 

• A null apparent horizon at r = r_,t = oo, where r_ is the smaller positive root of 
/(t = oo) = 1 - 2m/r - H^r^ = 0. In the case i^o > we will demonstrate that this is 
a regular black hole event horizon. In the case i^o = it is a null singularity. In both 
cases it is a finite distance from points in the interior. 

• A null surface ending at the point r = r^,t = oo, where r+ is the larger positive root 
of f{t = oo) = 0. For Ho > this is a cosmological event horizon. For Hq = it moves 
off to r = oo, t = oo and becomes a null FRW infinity. 

• A spacelike apparent horizon at the smaller positive root / = with t finite. It is 
connected to . . . 

• ... an apparent horizon at the larger positive root of / = with t finite. The two 
branches of apparent horizon at finite time t link up at a bifurcation point and delimit 
a normal region of spacetime from an antitrapped region. In asymptotically de Sitter 
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space, the outer branch of the apparent horizon asymptotes to the null surface through 
r = r+,t — !■ oo, whereas in asymptotically decelerating FRW cosmologies it ends at 
timelike infinity [15]. While the coordinate t is everywhere timelike, the coordinate r 
is spacelike only inside the normal region, enclosed by the apparent horizons. 

The casual structure of the McVittie geometry is shown in Figure 2. 



Figure 2: The conformal diagram of the McVittie solutions which asymptote the 
Schwarzscschild-de Sitter geometry in the future. The red dashed line is the surface 
r = r_,t = oo, and the blue dashed line is the cosmological horizon through r = r+. 
The broken curve at the bottom is the inhomogeneous Big Bang at r = 2m and finite times 
t, and the thin green dashed line is the union of the two branches of apparent horizons for 
finite t. Geodesic expansion is indicated by the now-standard 'Bousso wedges'. 



4.1 McVittie's Apparent Horizons 

A feature of the geometry often helpful in understanding its causal structure is the apparent 
horizon, a surface where at least one congruence of null geodesies changes its focusing prop- 
erties (see, e.g [2, 16]). As it crosses the apparent horizon this family of geodesies flips from 
converging to diverging (or vice versa). So to find the apparent horizon, we seek the locus 
of a vanishing geodesic expansion. For any observer, this surface is the light front on which 
the area is extremal. In spherically symmetric spacetimes, the light fronts are spherical, and 
so it suffices to look for the extrema of the spherical area around the observer along radial 
null geodesies. These satisfy the null constraint ds"^ = in (7), which yields 



where the upper (lower) signs apply for the outgoing (ingoing) rays. 

The spherical area is just A = 47rr^, so its variation (8) is SA = SirrSr = IGnrj^^SX, 
where A is the affine parameter along the null congruence, and ^ is the inverse of (8). Since 
we are interested in all the apparent horizons, for arbitrarily oriented null geodesies (i.e. 




dt 
dr 




±1 



(8) 
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both in- and out-going), this equation shows that to find it we should find the composite 
hypersurface for which A' = either along ingoing or outgoing geodesies. Thus, the full 
apparent horizon is defined as the hypersurface on which 




0, (9) 



which yields 



/ = i-^-i/(t)V = o. (10) 

To find out how each individual family of geodesies behaves on and beyond the apparent 
horizon we need to be a bit more formal [16]. Define the geodesic expansion by considering 
the Lie derivatives of the area two-form perpendicular to the geodesic congurcnce. The 
two geodesic congruences are given by the fiows of the vector fields dt + (^)_|_ <9r, the area 

form for spherically symmetric wavefronts is just sin 6d6d(f), and so the Lie derivatives are 
C±ui — 9 ±00, where the geodesic expansions 9± are 




(11) 



Where 9±^s are positive, the corresponding family of geodesies are diverging, and conversely, 
where they are negative the geodesies are converging. The spacctime is naturally divided 
between the regions where both expansions are positive (anti-trappcd) or negative (trapped), 
or they have different sign (normal). The question is, which of these regions may be accessed 
by observes who, at some point, can communicate with each other, and not have to face the 
doom of singularities separating them. 

In the limit of the Schwarzschild-de Sitter geometry, H — const, the apparent horizon 
coincides with the event horizon (as is clear from (10)). For a general McVittie solution the 
apparent horizon is given by the roots of (10), which express r as a function of t through 
the time-dependence of H. Evidently the area of the apparent horizon at any given time t 
coincides with the Schwarzschild-de Sitter horizon area for a spacetime with the curvature 
given by the Hubble length at that t. As time goes on the area of the apparent horizon 
changes, in such a way so that the smaller root of (10) moves inward as H decreases with 
time towards its asymptotic de Sitter value, while the larger root moves outward. 

We see this quantitatively from taking the time derivative of (10) and evaluating the rate 
of change of the radial extent of the apparent horizon. It is given by 



dr 
di 



2HHr^ 



horizon {l-V3Hr){l + V3Hr 



(12) 



where we have used (10) to eliminate terms oc m/r. Now, r = is the value of the 
horizon radius for the extremal Schwarzschild-de Sitter solution. So we see that at small t 
when H > the apparent horizon does not exist. As t increases its two spacclikc branches 
bifurcate and, as noted above, move inward (to smaller r) and outward respectively. The 
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inner branch evolves along spacelike directions, since the interval between its subsequent 
locations is 

2m 

using 1 = H^r^ on A.H. , (13) 

r 

and noting, from (12), that f < 0, such that (is^ > 0. The outer branch eventually becomes 
timelikc, asymptoting to null in the distant future. 

So far we have investigated the solutions r{t) oi f — with t finite. But in the limit 
i — > oo at least one of the two solutions remains finite. Let r± be the two positive roots of 

/(^ = oo) = 1 - ^ - Hy, (14) 

where Hq is the asymptotic future value of the McVittic Hubble parameter. The solution 
r — r-,t — oo is evidently a null apparent horizon, since in this limit f — 0, H — 0, and 
the norm of the normal vector to the surface vanishes. The properties of this null surface 
are key to our analysis. The pressure singularity in (4) may be absent there, since H = 0. 
To check, we must determine the rate at which H asymptotes to its constant future value 
along infalling geodesies, and at the same time settle the question of whether ingoing null 
geodesies reach it at finite affine parameter. We will do both in the next section. 

The larger root r = T+,t = cxo is a point on the cosmological event horizon, defined by 
the past-directed, outgoing radial null geodesies that end at it, when Hq ^ 0. Near it the 
geometry is controlled by the positive cosmological constant, and future-directed outgoing 
null geodesies are well behaved, crossing the horizon at finite values of the cosmic time 
t. Constant t surfaces also cross it (since they are spacehke). In the hmit Hq — 0, this 
horizon goes to infinity and becomes FRW null infinity. Note, that the finite-i branch of the 
apparent horizon at finite times — given by the larger root of (10) — always remains inside 
this cosmological Schwarzschild-de Sitter horizon, asymptoting to it at late times. Thus, 
the apparent horizon becomes first timelike and then eventually null, as noted above and 
sketched in Figure 2. 

To complete our description of the the causal boundaries of the McVittie solutions, 
we need to see what happens in the past. In the region covered by our coordinates the 
cosmologically expanding McVittie spacetime (1) starts from a past spacelike Big Bang 
singularity at r = 2m and approaches the de Sitter limit at t — )■ cxd if i7 ^ const 7^ 0, 
or some power law FRW cosmology if if — > 0. As already noted, the r = 2m, t — finite 
singularity is spacehke. It appears to the future of the naive FRW singularity at a = 
(where H diverges). 

One way to see this is to recall that r = 2m corresponds to a(t)|x| = m/2 in the 
coordinates (1), which yields t ~ l/|a;|^/'* for the singular surface at early times (where s > 
depends on the matter equation of state). Thus the singular surface at r = 2m approaches 
the surface a{t) = at spacelike infinity from above. Since in the de Sitter-energy fluid 
geometry the constant t surfaces asymptote towards the central observer's particle horizon 
as X ^ 00, the fi = 1 surface bends up in a similar fashion far away (see, e.g. [15]). Closer 
in, this surface also gradually bends up (relative to a{t) = surface), while remaining both 



ds^ = {dtf 



dr 
di 



-2 + 



1 dr 
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spacelike and below the apparent horizon at all finite times, eventually crossing the null 
surface r = r_, t = oo. 

So we see that the McVittic solution (1) with Hq > describes the region bounded 
by the null surface r = r_,t = oo, the big bang singularity r = 2m, and the asymptotic 
cosmological region. We can map it onto the Schwarzschild-de Sitter conformal diagram as 
illustrated in Fig. (2), where the jagged line is the big bang. Manifestly the spacehke big 
bang and the cosmological sector in the upper right are a part of the spacetime, since the 
null geodesies (8) get to the cosmological horizon r = r+ in finite affine parameter. 

In the next section we will concentrate on the null surface r = r_,t = cxo, first showing 
that it can be accessed in finite affine parameter by ingoing null geodesies, and determine 
whether it is traversable or singular. 

5 r = r_, t = oc: a Surface Too Far? 

In a typical picture of black hole formation from collapsing matter the apparent horizon is 
always found to be behind the event horizon, when matter satisfies the null energy condition. 
In contrast, in the McVittie geometries the apparent horizon extends outside the would-be 
event horizon. This might tempt one to argue that the expanding McVittie solution cannot 
have any black hole event horizons, since the future-oriented ingoing null rays seem to cross 
the inner branch of the apparent horizon. Led by this, Nolan [2] decided that the null surface 
r = r_,t = oo is a null boundary of the spacetime, separated from any interior point by an 
infinite lapse of affine parameter along any causal geodesic, and that therefore the McVittie 
spacetimes are geodesically complete to future directed null geodesies. 

Such situations are sometimes encountered as, for example, in extremal magnetically 
charged black hole spacetimes [17, 18] and some charged 3-D black strings [19]. However, 
they tend to be 'unstable' in the sense that when another parameter is turned on, such as 
e.g. electric charge, the boundary moves inward and the spacetime becomes geodesically 
incomplete. Hence it would be surprising if the surface r = r_,t = oo is a boundary in 
a general case, since this conclusion would have to hold for any external matter equation 
of state, in a completely generic fashion. Already in [2] doubts were expressed about this, 
although apparently the author decided against them. 

The argument of [2] can be invalidated quickly once one recalls that the ingoing geodesies 
that start at the Big Bang singularity and cross the apparent horizon in fact enter the normal 
region of spacetime. They are converging. The outgoing geodesies are not, but this does not 
say anything about the surface r = r-,t — oo. In fact the formulas for geodesic expansion 
(10) show that it is actually a null branch of the apparent horizon, and there a family of 
geodesies will change its convergence properties. It may happen that both are convergent 
beyond this surface, which would indicate a trapped region, and identify r — r-,t — oo as 
the event horizon. Thus to investigate the nature of the surface one must check explicitly if 
null ingoing geodesies can access it in finite affine parameter. 

In what follows we shall prove explicitly that they can, and hence that the surface r = 
r_, t = oo is not a boundary and the spacetime is not geodesically complete to future ingoing 
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null rays. We will not be able to explicitly solve the geodesic equations.^ However an exact 
solution is not necessary because we can set a finite upper bound on the lapse of the affine 

parameter along ingoing null geodesies during the approach to r = r_,t = oo. 

To proceed wc need the null radial geodesic equations. One is the first integral (8). 
We get another by varying the geodesic "Lagrangian" C = . !1 , "j^"^^ * — ft'^, where 

i zm/r a/1 — 2m/r 

prime denotes derivative with respect to the affine parameter. If we define two functions 
F± = ±(yT^^^rrzyV =F Hr)t' + r' / y^l — 2m/r, the Lagrangian factorizes as C — F+F_. This 
means that the geodesic equations will be of the form 

(F+a,.F_ + F_dq,F+)' = F^d,F_ + F_d,F+ , (15) 

where q is either t or r. Specializing to ingoing null geodesies implies by (8) that we can 
set F^ = after the variation, reducing (15) to only the terms involving derivatives of F_. 
Choosing q = t and using (8) gives 

dt _ 1 

dr ^1 - 2m/r{^l - 2m/r - Hr) ' 

r" = — (16) 

^1 - 2m/r(^l - 2m/r - Hrf ' 

where the oc t' terms in the second equation are eliminated using the first. 

Prom Eqs. (16) we can see by a simple argument that ingoing null geodesies reach the 
surface r = r_, t = oo at a finite value of the affine parameter. Combining Eqs. (3) and (4) 
gives 

\^^H = -^t:Gn (p + p) < , (17) 

where the inequality holds if the null energy condition is satisfied. In this case the geodesic 
'acceleration' r" in the region ji <1 (the future of the big bang) is non-positive, and therefore 
an ingoing geodesic with initial geodesic 'velocity' Tq < can never turn around or decrease 
its speed. Starting at any initial distance r > r_, it must reach r_ after an affine 'time' 
AA < (r — r_)/|rQ|, which is finite. The first of Eqs. (16) shows that in this limit t — ^ oo 
as expected. This conclusion is completely independent of the asymptotic form of H as 
long as the null energy condition is satisfied. Therefore, all physical McVittie solutions (1) 
are geodesically incomplete to null future oriented ingoing geodesies, in contradiction to the 
claim of [2]. 

This suffices to show that the McVittie spacetime (1) is geodesically incomplete, but it 
can be refined in order to produce a more precise estimate of the rate at which the surface 
r = r_,t = oo is approached. We will need this estimate in the analysis of when the null 
surface r = r_,t = oo is traversable and is a horizon, and conversely when it is impassable 
and a singularity. Wc defer the details of this analysis to the Appendix A. 



^For example, in the case of single fluid solutions and vanishing cosmological constant, where a = V, the 
geodesic equations reduce to a single first order Abel's differential equation ^ = ^Z+fS)^^"'" ^tT^^ 
with / = which, as far as we could ascertain, is not known how to integrate [20]. 
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6 In Medias Res: the McVittie Event Horizon 



So, what is the null surface r = r^,t = oo? We will show that in the case when H — )■ 
Hq > 0, it is a regular event horizon, and the McVittie solution really is a black hole, which 
asymptotes the Schwarzschild-de Sitter solution. To prove this, ideally we'd construct a 
family of ingoing null geodesies that cross the horizon. This family is parameterized by an 
integration constant u, and each line depends on the affine parameter v = X. Given such a 
construction we could view the smooth parametric equations of the worldines 

r — r{u,v), t — t{u,v), (18) 

as a coordinate transformation, exchanging the original coordinates r, t (with t divergent and 
thus ill-defined on the horizon) with the new coordinates u, v which arc perfectly regular and 
well defined on the horizon and beyond. We cannot find such coordinates globally in closed 
form due to the complexity of the nonlinear geodesic equations. However, we will now 
demonstrate that such coordinates do exist, and we can construct them at least locally in 
the near horizon region. Further, once we have demonstrated that those coordinates exist, we 
will check that near the horizon they reduce to null coordinates which cover Schwarzschild- 
de Sitter, and use this to prove that all curvature invariants are finite at the surface, in 
contrast to McVittie solutions which have asymptotically vanishing Hubble parameter.^ We 
will demonstrate this at the end of this section. 

First off, a set of null coordinates which cover (7) certainly exists, since by spherical 
symmetry the r — t part of the metric is conformally fiat. Next, the metric (7) is clearly 
non-degenerate at r — r^,t — oo when Hq > 0, since at / ^ it reduces to 

ds'^ ^ -2drdt + + r^_dn2 . (19) 

At this surface the coordinates r,t are 'almost' null, in the sense that the only term which 
spoils the null tilt of the t coordinate is dr"^ / Hlr"^, and in addition there is the irritating fact 
that t diverges. However, the individual curvature components remain finite there. ^ 

The unfortunate features of the coordinate t can be corrected with a radially dependent 
time translation. Choose a new time coordinate 

dv 

dT^dt+ ^ , (20) 

Jl — 2m/r(Jl — 2m/r — F{r)r) 

where F{r) = H{t{r)), and t{r) is a one-parameter family of solutions of Eqs. (16) describing 
ingoing null radial geodesies. This new time coordinate is well defined away from the extremal 

^Checking that the curvature invariants are finite on a nuU surface is not enough to conclude that the 
surface is regular. An example is provided by wavy cosmic strings with null singularities which are invisible 
to any invariants [21]. However, our proof is stronger since wc demonstrate the existence of a smooth family 
of ingoing null geodesies. The calculation of the curvature invariants then simply serves to show that there 
are no hairs on this surface, and that the geometry really asymptotes Schwarzschild-de Sitter there. 

^Curvature components arc easy to calculate using either (1) or (7). We give them for the r, t coordinates 
and the metric (7) in Appendix B. One can verify that they are regular at r = r_,f = oo, and by the 
'staring' argument one can begin to suspect that any number of the covariant derivatives of this metric will 
remain regular... at least as long as ^ i/o 7^ on this surface. 
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limit, and can be integrated directly, because the second term on the RHS of (20) is a total 
derivative. Also a suitable coordinate exists in the extremal Hmit as well. This new time 
coordinate does not diverge near the surface r = r^,t = oo! Indeed, in terms of this 

coordinate the first of Eqs. (16) reduces to 4^ = , ^ ( , ^ , ^ ) , 

^ ^ dr ^l-2m/r \ ■s/l-2m/r-F{r)r ■>J\-1mlT-H{t)r ) ' 

and so along an ingoing null geodesic, where we insert the function tif) in place of 

At = 0. (21) 

Thus r is a conserved quantity of a given null geodesic, and can only vary from one geodesic 
to another. Now, note that 

• the transformation of coordinates (20) is nonsingular everywhere away from r = r_; 

• the ingoing null geodesies reach r = r_ in finite affine geodesic parameter, as we have 
shown in the previous section and in more detail in the Appendix A; 

• the family of ingoing null geodesies is smoothly parameterized by the affine parameter 
and an integration constant which indexes the geodesies. 

This implies that the coordinate transformation (20) remains nonsingular in the whole 
patch of the exterior geometry inclusive of the surface r = r_. Because of this, the coordinate 
T, being orthogonal to the affine parameter A since it is conserved along any given geodesic, 
can only be a function of the integration constant in the solution (16). This implies that r 
is continuous, and so also finite, as claimed. Thus, the divergence of the time coordinate t 
has been regulated away, and r will be well-behaved - smooth and bounded - on the surface 
r = r_, t = oo. 

We can now rewrite (7) in terms of the r coordinate. After some straightforward algebra 
we obtain 

ds- = -fdr-- [et- , ^ 

^1 - 2m/r \ - 2m/r - F(r)r^ 

|+.^«„(22) 



l-2m/r\^ l\ _ 2m/r - F[r)r U\ - 2m/r - F{r)rf 



where / = 1 — 2m /r — H^r^ is still the same function / we defined in (10), and H now 
depends on both r and r, through the integral t{r,T) defined by (20). When we take the 
limit r = r_,t = oo, thanks to the form of F{r) — H{t{r)) and the limiting form of the 
function / when r — > r_, this metric reduces simply to 

ds'^ = 2dTdr + r^dOs , (23) 

which is non-degenerate. The coordinates r, r are regular and bounded at this surface. 
Further, r is defined by a solution of differential equation (20) which exists, and is analytic 
everywhere in the domain of interest. So the full metric of Eq. (7) must also be analytic at 
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r = r_, i = oo, when H ^ Hq 0, and the surface r — r^,t — oo must be a traversable 
horizon. The geometry (1) really represents a black hole embedded in an expanding universe, 
which emerged from some singular, inhomogeneous initial data at the Big Bang singularity 
= 1, so long as Hq > 0. 

What goes wrong with this argument when i?o = 0, as for example when a — P'l The 
surface r — r_ — 2m is still null and still at finite affine distance. However a glance at Eqs. 
(20) and (22) shows that the coordinate r is not finite on the "horizon" in that case. As we 
will see, that is not due to a bad choice of coordinates — the surface is actually singular. 

To see this we begin by examining the covariant curvature components, Christoffel sym- 
bols, and inverse metric given in Appendix B. One can see immediately that the only 
singularities in these components occur when r = 2m and when a{t) — 0. Taking (ordinary) 
partial derivatives with respect to i or r can change the degree of divergence at these points, 
but cannot introduce divergences at any other points. Since all possible curvature invari- 
ants are constructed by products of these functions and their ordinary partial derivatives, 
divergences can only occur on these two surfaces. 

We have already seen that the surface r = 2m, t — finite is a big bang curvature singu- 
larity. The surface t = (where a vanishes) is to its past and not a part of our coordinate 
patch. What about the surface r = r_,t = oo? Since r_ > 2m when Hq > 0, no curvature 
can diverge on the horizon, and therefore it is traversable and regular. 

To evaluate the curvature on this surface in more detail we should take the limit as we 
approach it along an ingoing geodesic. For the case Hq > 0, the near horizon scahng is 
t{r) ~ ^In^:^^, where a > 0, as we have shown in Appendix A. Looking again at the 

ingredients of the curvature invariants, recalling that H ~ e~^°* and r_ > 2m, one can see 
immediately that the only terms that survive in the limit are those proportional to H and 
not to any of its time derivatives. But H = Hq on the horizon, and therefore all curvatures 
on the horizon are exactly equal to their values on the horizon of a Schwarzschild-de Sitter 
black hole of mass m and Hubble constant Hq. 

In fact, the observation that the inverse metric, Christoffel connection, the curvature 
tensors and all of their derivatives are regular on this surface in effect guarantees the existence 
of a nonsingular family of geodesies which cross this surface, since one can start in any 
arbitrarily small neighborhood of the surface and construct the geodesies using the power 
series method. This is why we were able to define the coordinate r in (20) and change the 
metric to (22), which is regular across the horizon. 

The case Hq = 0, a{t) — )> t^ at large t is qualitatively different. Since r_ = 2m, the above 
argument no longer guarantees that the curvatures arc bounded on the would-be horizon. 
On the contrary, we find that they are in general divergent. The equation for an ingoing 



null geodesic in the near-horizon region (with a{t) ~ at large t) is t ^ pr/^Jl — 2m/r. To 
determine this scahng, susbstitute a f^H — p/t\n the first of Eq. (16). This gives 



Because the null geodesies do reach the surface r = r_,t = oo in finite affine parameter, t 
must diverge to infinity as r — > r_ = 2m when i^o = 0, as above. Indeed, if this were not true. 




dt 
dr 




1 



(24) 
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and t remained finite as r — > 2m, null future oriented ingoing radial geodesies would have all 
ended on the Big Bang singularity - a past spacelike singularity - which is clearly impossible. 
Further, the sign of ^ cannot change as t increases, because if it did it would imply that 
the future evolution would cease, with null geodesies turning back in time as r ^ 2m. This 
would imply the existence of closed timelike curves and breakdown of causality in the regions 
of McVittie solution (1) which are obviously fine, and so it cannot happen. Hence, we see 

that 1/t must follow ^^'^^^^ from below. Setting e = ^J^^^-2mJr — pr/t > 0, solving for 
t and taking its derivative to the leading order in e, and comparing it to the formula (24), 

that translates to ^ = — (e-^l — 2m/r)~^, we can exphcitly confirm this scaling, when we 
pick 

e = A(1- 2m/r) + . . . , (25) 

for A a constant and ellipsis designating subleading terms. 

Using this scaling one sees that the Ricci scalar R = 12H^ + 6H / — 2m/r is finite 
in the limit. However one may suspect that higher curvature invariants will not be, since 
by taking a sufficient number of partial derivatives with respect to r one should be able 
to obtain a divergence from the square root term. To be certain there are no miraculous 
cancellations, we computed the invariant S = (V ^VuRpaSx)'^ (see Appendix C for the result). 
It contains terms such as /J*^(i/)^(1 — 2m/r)~^ which diverge at the horizon along ingoing 
null geodesies, and there is no cancellation among the divergent terms. 

An alert reader may by now be a little discomforted by the following observation. On 
the one hand we are claiming that all curvatures are finite on the horizon when Hq > 0, and 
have a finite limit when Hq — > (since the horizon curvatures of a de Sitter-Schwarzschild 
black hole of mass m simply reduce to those of a Schwarzschild black hole of the same mass 
in that limit, which are themselves manifestly finite). On the other hand we claim that 
certain curvature invariants diverge when Hq — 0. How can this be? 

The resolution of this puzzle is that some curvature invariants, such as S, do not mono- 
tonically decrease along ingoing null geodesies when Hq > 0. Instead, S reaches a maximum 
a finite distance from the horizon r = r_ and then decreases to its Schwarzschild-de Sitter 
value at r = r_. As Hq approaches zero, this maximum increases and at the same time its 
distance from the horizon decreases. In this manner the curvatures on the actual horizon 
r — r- are bounded (and small if m is large) for any Hq > 0, but nevertheless the horizon 
becomes a null singularity at Hq = 0.^ The order of limits where the cosmological constant 
is held fixed and one goes to the horizon at t = oo, and then sends Hq — ?> 0, regulates 
the singularity by first "infiating away" the FRW matter and pressure to zero due to the 
cosmological expansion.^ No such mechanism is in operation when Hq = from the start. 
Nevertheless as we have seen explicitly, the singularity that is present in the case oi Hq — 
is very soft: first of all it is null, and secondly it is weak, being relegated to invariants which 
involve at least two derivatives of the Riemann tensor. Strictly speaking it does not violate 
the spirit of the strong cosmic censorship conjecture (see e.g. [22] for a nice review). 

'^This is reminscent of the Gibbs phenomenon encountered when one studies discontinuous Umits of Fourier 
series representations. 

^Old age often leads to hair loss. 
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7 Summary 



In this paper we analyzed the McVittie solutions of Einstein's field equations, which have 
been debated for nearly 80 years as candidates for describing the gravitational fields of 
spherically symmetric mass distributions in expanding FRW universes. For simplicity we 
focused on spatially flat McVittie geometries. Our results show that the McVittie solutions 
that asymptote to FRW universes dominated by a positive cosmological constant at late times 
are black holes with regular event horizons. While we were not able to construct the explicit 
form of the extension of the geometry past the black hole event horizon r = r_,t = oo, we 
showed that such analytic extensions exist when H ^ Hq > 0, being defined by the ingoing 
null geodesic family. Further, we demonstrated that in this limit the near horizon geometry 
is asymptotically Schwarzschild-dc Sitter, which is essentially a consequence of the cosmic 
and black hole no-hair dynamics. Thus at least in this case, the shaded area past the surface 
r = r-,t = oo in Figure 2 is a real black hole interior of some sort. 

When the exterior geometry asymptotes to an FRW cosmology with a power law scale 
factor, so that H — p/t ^ in the future, the McVittie solutions are singular on the 
null surface r = r_ = 2m, t = oo. This surface can still be reached by null geodesies 
in finite affine parameter, and so is a real singularity. However it is very soft since it takes 
invariants involving at least two derivatives of the curvature to detect it. While it is naked, it 
respects the spirit of the strong cosmic censorship [22] in that it cannot be observed directly 
by external observers. Our results correct the claims made in [2] which have propagated 
through the literature.^ 

Our most interesting conclusion is that the McVittie solutions in spacetimes which are 
asymptotically dominated by a positive cosmological constant are black holes. It seems rather 
remarkable that these black holes were not identified until nearly 80 years after the metric 
was written down. One might be puzzled that what seems to be an infrared quantity — the 
cosmological constant — regulates a short-distance property of the geometry, the would-be 
null singularity. This can be understood as a consequence of the special and restrictive form 
of the McVittie solutions, which requires the validity of a continuous fiuid description down 
to the shortest scales, prevents the hole from accreting, and enforces that its energy density 
be exactly homogeneous on some set of spatial slices. As we have discussed, these conditions 
are not particularly well-motivated physically. One could give up spatial homogeneity of the 
energy density, as in e.g. [12] and its generalizations, or allow the hole to accrete. 

Another strategy for describing black holes in cosmological backgrounds is to go beyond 
the simple fluid description of cosmological matter and replace it with a microscopic descrip- 
tion. Such ideas motivated the work of [9], and were behind the proposal of the levitating 
dark matter in [24]. A concrete example which realizes this has been found very recently 
in [25], based on the papers [26, 27], and generalized in [28]. The solutions found in those 
papers are very interesting, in that they encode large distance properties of the McVittie 

^See for example [23]. There are many conflicting incorrect claims, ranging from assertions that McVittie 
solutions are geodesically complete, over the suggestions that they cannot be black holes because of exotic 
spacelike singularities, to the statement that McVittie solutions in any accelerating universes are black holes. 
Our results establish a systematic method for determining what the McVittie solutions are and dispel those 
incorrect claims. 
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family. The matter content is a scalar field with an asymptotic runaway potential and a 
system of gauge fields with nonzero charges. Far from the hole the scalar rolling down the 
potential provides the source that drives the cosmic expansion. Near the hole the charge 
contributions correct the effective potential for the scalar and give it a large mass, employing 
essentially the same dynamics as the supersymmetric attractor mechanism in asymptotically 
flat black holes [29]. Thus close to the black hole horizon the scalar decouples and there 
is no medium that could create divergences. This method of regulating the null McVittie 
singularity is local, but of course it requires abandoning both the fluid description near the 
black hole and the insistence of homogeneity of energy density outside of it, and is therefore 
completely consistent with our findings. It suggests that there might be other mechanisms 
for regulating the near horizon geometries of solutions which reduce to McVittie far away 
that might warrant further study. 
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A The Short Happy Fhght to r = r_, t = (X) 

Here we give the refined argument that the ingoing null geodesies reach the surface r = 
r-,t = oo in finite affine parameter, which makes the McVittie geometry of (1) geodesically 
incomplete to them. We further compute the precise rate at which the null radial geodesies 
approach the surface r = r_,t = oo, which we will need to show that in the limiting case 
with if 7^ this surface is a regular event horizon. 

So to start with, note that the first of the Eqs. (16) has a pole on the right-hand side 
there, which we can see by rewritting it as 

— 9m I r -\- I-fr] 

(26) 



dt r (^^1 - 2m/r + i/r) 

dr H'^^Jl — 2m/r{r — r+)(r — r_)(r + r+ + r_) 



where r± are the two positive roots of (10), and — (r+ + r_) is its one negative root. In the 
limit r — >■ r_, this yields 

dt — >■ (27) 

where we in principle must take into account the fact that the roots r± of Eq. (10) are 
time-dependent, approaching the roots r± of (14) only as t ^ oo. 
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Equation (27) defines a differential equation that controls how the geometry approaches 
the surface r = r_,t — oo. Since we are interested in the limit r — )■ r_, where the only 
divergence in (27) appears as the denominator of the RHS vanishes, while all other quantities 
remain finite, we can set r± to their asymptotic values everywhere except in the difference 
r — r_(i). This yields a first order linear differential equation 

'i-^{r-r-i*)), (28) 

where from (27), we find a = ^){r++2r ) • -^^^J from the extremal limit where r± 

coincide, and when the McVittie geometry asymptotes to the Schwarzschild-de Sitter limit 
with H — >■ Ho ^ 0, integrating this is straightforward. To the leading order, when m <^ H 
(with appropriate generalizations as m approaches the extremal value), in the limit t — )■ cxd 
we have r_(t) — )■ r_ + Ae~*'^"*, where A = 8m"^p(t = 0). This follows from Eq. (3) in 
the limit t — )■ oo and m <^ H, substituting for the energy density contents of the universe 
p = A + po (^) , where s — 3(1 + w) and w — p/p is the energy fluid equation of state, so 
that H ^ Ho + 0{e-'"°*). Therefore the solution of (28) is 

(1 - sa)r_e~^"*/" + Ae"^^"* ^ (1 - sa) (r - r_) + . . . , (29) 

to the leading order, where the ellipsis denote subleading terms, and a > reflecting the 
fact that t ^ oo as r ^ r_. Clearly the rate at which r approaches r_ depends on the cos- 
mological contents, since it is controlled by 1 — as. When as <1 the evolution is dominated 
by the second term on the LHS, while when as > 1, the flrst term takes over. Further, 
when as — 1, the LHS of the solution (29) degenerates into a resonance-like polynomial 
times the exponential e~^"*, and while this case needs to be given special treatment, its 
generic behavior is the same as the nondegenerate cases. Moreover, one can check that in 
the extremal limit r+ = r_, which occurs when m is very large, the right-hand side of (8) 
has a double pole, implying that the divergence of t is even faster, going as 

Hot ^ ^TTT^ r . (30) 

3Ho{r - r_) 

By the arguments above, H oc e~*^°* in this limit. Note, that H < 0. Hence, along null 
ingoing radial geodesies near the null surface r^,t — oo, we have if oc (r — r_) when as < 1 
(with essentially irrelevant logs when as = 1, which we shall ignore) and if oc (r — r_)"^ 
when as > 1, and thus the coefficient of r'^ on the right hand side of the second of Eqs. (16) 
behaves as 

^ y-C\r-r_y, (31) 

^1 - 2m/r(^l - 2m/r - Hr^ 

where p = —1 when as < 1 and as — 2 when as > 1, and is a positive constant given 
by the ratio of the energy fiuid density and the cosmological constant. This expression 
is generically divergent for typical values of as. The asymptotic form (31) enables us to 
integrate the second of Eqs (16) along a null geodesic, to find that to the leading order 

r' = r: exp ( -C^ f dr l'—!:^^ \ , (32) 
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where r'^ is the initial radial 'velocity' of an ingoing geodesies starting at some initial distance 
r = r* > r_. So wc confirm that the geodesic 'velocity' r' is non-positive, and it might vanish 
when r — >■ r_. Integrating the exponent in (32), we find that when as < 1, so that p — —1, 



'=<{v^r-' (33) 



which clearly implies that r' diverges to — oo as r — > r_, since r'^ < 0. Therefore, there is no 
stopping of the fall towards r_. Indeed, since the affine parameter is A = ^, we find that 



7 / dr( ) , 34 

1 Jr. ^r* — 



and since the integral is convergent, it takes a finite affine time AA to get to r = r_. Similarly, 
when as > 1, so that p = as — 2, we have 

^'-^exp(-^(^-^)""'), (35) 

where A = exp (^-^^^zi ( j*r^~ ) So ^ is finite and strictly negative, since r'^ < 0. 

Therefore, since A — * , we find that 



^AA 



Since the exponent vanishes as r — > r_, the integral still converges, and again it takes a 
finite affine time to get to r = r_. Hence, in any case the right-hand side of either (34) 
or (36) is always finite, confirming our assertion that the average geodesic 'velocity' r' is 
negative definite. A similar proof exists when the asymptotic value of H vanishes. This 
detailed analysis confirms our earlier argument, showing that ingoing radial null geodesies 
reach the null surface r = r_ , t = oo in a finite lapse of the affine parameter, and hence that 
the portion of the McVittie geometry (1) bounded by this null surface, the singularity /j, — 1 
and the exterior cosmological regime is not geodesically complete. 



B The Curvature and Einstein Tensors 

It is straightforward to compute the full Riemann tensor for the McVittie solution (1). 
Specializing to the coordinates r,t, and starting with the metric (7), after a straightforward 
albeit tedious calculation we find that the nonvanishing components are 

R - ^ 

^trtr I Q ) 

^Jl - 2m/r 

Rtote = - 2m/r + ^(l - ^) - H'r'{l - ^ - H'r') , R,^^ = R,ere sin' 9 , 
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R. 



_ H'^r'^ - m/r _ . ^ 

rdrO — —. ^ , Jrir<l>r<l> — J^rdrO Sm d . 

1 — 2m/r 



(37) 



/2m tt2 2\ 2 ■ 2 n r> Hv {H'^v'^ — m/ v) .2/, 

Rocbed) — [ \- H r )r sm , RtOre — / — , Rtsrs — RtOre sm , 

^ r ' ^1 - 2m/r 

and it is quite clear upon direct inspection that they are all well behaved in the limit 
T — T-^t — oo d& long as if — > ifo 7^ 0. 

The Einstein tensor is found from the contractions of the Riemann tensor. It is 



--H . 



(38) 



C The Singular Invariant in the Power-law McVittie 

Here we quote the explicit form of the invariant 5 = {S/ uRpa8\f"- 



2m 



12y 1 - —-r'°H^(-l?,2Qm^{2m - rf + r^(885m^ - 1686mV + 1240m'^r^ - 412mr^ 
+52r^)if2) - 24r^^ (-158m^ + 185mV - 63mV^ + mr^ + 2r^) H^HH 
-24(2m - ry°HH(m^ (6m^ - 19mr + 8r^) H + 



2m 



' 1 - — (67m^ - 56mr + 12r^) 
(l4m^ - llmr + 2r^) + 4r^(-2m + r)ii"^ -6mr^° (57m^ - 47mr + lOr^) 



-2mVl - (334m^ - 347mr + 96r^) + 3J1- ^(240m^(25m - 12r)(-2m + r)' 



+r" (lOW - 88mr + 19r^) ff^) - 6mr^if (l80m(-2m + r)^ + y 1 - —r\-5m + r)d^H 
+4(-2m + r) 2mV^° (6m^ + mr - 2r^) + 3^1 - ^r^^ (57m^ - 52mr + 12r2) 
+r^(-2m + r)ii'2 ^771^! - ^ (847m^ - 832mr + 222r2) - 6(5m - 2r)r^9i''i7 j 

-2m2(2m - 1 180m(2m - r f + 7 Jl - —r^d^H 



+1I 



2m 



-2m + r) (-40mV^^if2 + 3(-120m2(2m - r)^ (oSm^ - 60mr + Mr^) 
(ri^(l - 2m/r)ii/^)-i . 
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